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Abstract: The arrival time coherence of particles in the Ultra High Energy Air Showers 
where the center of mass energy of the interaction is of the order of IQi^^eV , puts strict 
constraint on the propagation of particles in a hypothetical extra-dimension. We first argue 
that at such high energies bulk modes and massive KK-modes can be produced abundantly 
and in many models their phase space volume is larger than confined modes. Then, we study 
the minimum propagation time in one and two-brane models and show that a large part of 
the parameter space of these models are ruled out unless the confinement of fields is proteced 
by symmetries up to energies not accessible even to the high energy tail of Ultar High Energy 
Cosmic Rays (UHECRs). As a by-product we confirm the result obtained in some previous 
works about the close relation between a small Cosmological Constant and the hierarchy 
problem. 
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1. Introduction 

Ever since the proposition by Th. Kaluza and O. Klein in 1920s to use a 5-dimensional 
space-time for unification of Gravity with Electromagnetism [1] , space-times with more than 4 
dimensions have been the hope of physicists to solve problems of High Energy Particle Physics. 
Last ideas in these series are suggestions by N. Arkani-Hamed, et al. [2] and by L. Randall 
and R. Sundrum [3] for using large extra dimensions and localized matter to solve mass 
hierarchy problem inspired by some previous works of V. Rubakov and M. Shaposhinkov [4] 
on domain walls in higher dimensional spaces and P. Horava and E. Witten [5] on M-theory 
models with Compactification in spaces with D-brane boundaries. 

In the first proposals only gravity could propagate in the bulk. It has been however 
found that the total localization of all fields except graviton on 3-branes is not realistic. In 
fact brane solution are cosmologically unstable and at least one scalar bulk field (radion) [6] [7] 
is necessary to stabilize the distance between branes. In some brane models inflaton [8] also 
has to propagate to the bulk to make inflation with necessary properties. A deeper insight 
to the propagation of gravitational waves and massive particles with bulk modes in models 
with infinite bulk has illustrated that even the warping of the bulk can not stop their escape 
from branes [10] [11]. 
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Most of localization mechanisms are evolutionary i.e. based on special configuration of 
matter fields with localized properties like topological defect solutions which can arise during 
phase transitions in the Early Universe [12]. In these models the real dimensionality of the 
space-time is larger than 4 but our Universe is confined to a domain wall (3-brane) where 
fields specially at low energies (with respect to quantum gravity scale) are geometrically or 
gravitationally localized. However, geometrical settings like a warped metric are not always 
enough to confine fields on the branes. It has been shown that in spaces with > 3 extra- 
dimensions gravity can not be geometrically confined to a 3-brane and p-form fields in the 
bulk must be added to stabilize the brane (defect) [13]. In 5-dim. models gravity and scalar 
fields can be localized on the brane with negative tension [14] (or the brane with smallest 
value of warp factor in 2-brane models of [15] [16]). Warp geometry can localize fermions 
only on the positive tension brane (which can not solve the hierarchy problem). Vector 
fields can not be geometrically confined. Localization of gauge fields and fermions on the 
negative tension brane is achievable through special particle physics setups [14] [17]. The 
localization scale however is considered to be not much higher than warping scale, otherwise 
a new hierarchy can appear [18]. At higher energies one expects that symmetry restoration 
(e.g. chiral symmetry of fermions) leads to escape of particles from the brane. 

Even when a warped geometry is enough to confine fields on the brane, the wave function 
of the zero mode can penetrate to the bulk (but has an exponential maximum on the brane) [3] . 
In infinite bulk models KK continuum begins from m = and this affects the long range 
behavior of gravity and other massless fields [10] [19] [20]. Massive fields if they have bulk 
modes can decay to the bulk with a life time which depends on the fundamental scale of 
gravity [19] [20]. For orbifoldized models the spectrum of KK modes is discrete. The long 
range effect of massive graviton modes is less important but the probability of decay of massive 
modes to the bulk is unchanged (see next section). Universal extra-dimension models in which 
all SM particles propagate to the compact dimensions are not ruled out for compactification 
radius of order TeV~^ or even lower [21]. 

At present brane world models have been constrained only based on the probability of 
direct observation of processes involving the production of gravitons and its Kaluza-Klcin 
modes [22] [23] [24]. A detail investigation of observable signal of the RS type models in 
Tevatron and LHC are performed in [25] and KK-mode production in the early Universe 
in [26]. The existent and near future accelerators can constrain the scale of gravity (and 
thus the size of the extra-dimensions) up to ~ 30Tel^. The interaction of Ultra High Energy 
Cosmic Rays (UHECR) with protons in the terrestrial atmosphere has a CM energy close to 
lOOOTey and is the most energetic interaction of elementary particles we can study today. It 
can be used to constrain the fundamental scale of gravity and the compactification scale up 
to much higher energies. 

The mass of KK modes detected by an observer on the brane is the result of smeared 
dimensions in the wave function. Classically however, it can be interpreted as a delay in 
the displacement of particles. For the observer on the brane if the delay slightly modifies 
the propagation of the particle in the detector, it is interpreted as a larger particle mass. 
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otherwise it is seen as an arrival delay, specially with respect to the particles which propagate 
only on the brane. In the case of an air shower, the time coherence of the showers will be 
destroyed. 

In this work we calculate the minimum propagation time of particles ejected to extra- 
dimensions for a number of warped brane models and compare them with arrival time res- 
olution of present Air Shower detectors. We restrict our study to D = 5 models. This is 
enough for understanding the general characteristics of the propagation in extra-dimension 
from the point of view of an observer on a (3 + l)-brane and can be considered as a special 
configuration for models with higher dimensions. Our attention is mostly concentrated on 
the classical structure of the brane models because results are independent of the detail of 
quantum field contents and origin of the branes. However, before doing this we must assess 
the possibility and the probability that UHECRs' interaction in the atmosphere can produce 
bulk modes i.e. escaping particles. Given the rarity of UHECRs, only models in which the 
probability of production of bulk modes is very high can be constrained by this method. 

2. Production of Bulk Modes 

The whole idea of constraining brane models with UHECRs depends on the possibility and 
probability that remnants of UHECR interaction in the atmosphere can penetrate into the 
extra-dimensions. In this section we review the localization of fields on the brane with a 
special attention on models which provide bulk modes i.e. not all fields are confined to the 
visible brane at all energies. We estimate the probability of producing these modes at energy 
scale of interaction between UHECRs and nucleons in the terrestrial atmosphere. 

By definition in universal models [21] any field has bulk modes and propagates in all 
space-time dimensions. The interesting case for solving the hierarchy problem is when the 
size of the compactified dimensions are of the order of weak interaction or lower^. This is 
~ 3 orders of magnitude smaller than the CM energy of UHECRs' interaction. Therefore in 
the case of universal models, UHECRs can produce lowest KK-modes abundantly. In [23] 
it is argued that UHECRs can not probe the physics at very high energies simply because 
their interactions is dominantly electromagnetic. It is true that probability of the exchange 
of a heavy particle e.g a massive boson related to symmetries beyond Standard Model is very 
small with respect to a low transverse momentum EM cascade. However, in universal models 
as SM fields have bulk modes at very high energies all dimensions are "seen" as to be the 
same and the parameter space of EM cascades with non-zero momentum component in the 
extra-dimensions is much larger than cascades restricted to the three infinite dimensions. In 
the language of KK-modes, with a good accuracy the lowest modes can be considered as 
massless and they can be produced with the same probability as zero modes. In the following 
subsections we argue that for non-universal models one expects that at the CM energy of 

^In some of universal models the extra-dimensions are not warped. Here we only study the propagation 
in warped spaces. Nevertheless, when the compactification scale is much larger than /x (see Sec. (^) for 
definition), the warp factor is very close to one and the results of following sections are applicable. 
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UHECRs interaction, most of localization mechanisms be no longer active and Standard 
Model particles can escape to the bulk. 

2.1 Scalars and Spin-2 Fields 

In warped models scalars and spin-2 fields can be confined geometrically. The zero mode of 
these fields has an exponentially decreasing wave function in the bulk [3] [27]. For models 
with infinite bulk the KK-spectrum begins from zero and therefore there is not a real confined 
zero mode. In addition, if the field has a non-zero 5-dim. mass, it has been shown [11] that 
4-dim mass eigen modes on the brane are complex and decay to the bulk with a width: 

r/m4 oc {rrii/^f ml = ml/2 (2.1) 

where m4 is the real part of the 4-dim. mass eigen value, m^ is the 5-dim. mass of the 
field, and /U is the warp scale in the static RS metric (Eq. ([4.1D below). If the fundamental 
scale of Quantum Gravity is comparable to the weak interaction scale, at CM energy of 
UHECRs interaction it is expected that due to radiative corrections even massless particles 
like gravitons have an effective non-zero mass. For short distances relevant to the propagation 
of UHECRs in the terrestrial atmosphere, massive modes have a Yukawa type potential and 
their coupling is exponentially suppressed [11]. However, if the m4 is much smaller than CM 
energy, the effect of exponential term in Yukawa potential is negligible. For most energetic 
UHECRs EcM ~ IQ^^eV. This means that the coupling to modes as massive as ITeV is 
roughly the same as massless modes. The width in (|2.lD depends on the effective 5-dim. mass 
of the field and the warping scale. We discuss their implication on the decay of massive modes 
to the bulk and on the test of brane models in Section ^. 

The above argument is also true in the case of 2-brane models where the spectrum is 
discrete. We show briefly that in static/quasi-static models the zero mode of scalar/spin-2 
fields with m^ > has an imaginary part i.e. it decays to the bulk (See also [28]). We 
determine the propagator of scalar/graviton using the Green function method discussed in 
detail in [27]. 

After changing variables y in metric ( [4.1D to: 

z = -e^'y (2.2) 

ds^ = dx^'dx'' - dz^). (2.3) 

we apply the boundary conditions to both branes. Without loss of generality we assume 
one of them is at w = or z = - = i? and the other at u = L or z = -e^^ = R' ^. The 
Green function (2-point propagator) A(x, z, j;', z') is the solution of 4-dim. mass eigenstate 

^We are only interested in the case where bulk is static. Therefore implicitly it is assumed that these points 
correspond to the fixed points of the radion field. 
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equation: 



{z'^dl + zd, + p^z'^ - d^)Apiz, z') = ^b{z - z') (2.4) 
A(x,z,:e',z') = j -^e'P(—')Ap{z,z') (2.5) 

^ (j) Ap(z,/) (2.6) 



d=^JA + R^^ml (2.7) 
The boundary and matching conditions for right and left propagators: 

A<^ = Ap{z,z') z < z' , Ay = Ap{z,z') z > z' (2.8) 



are as followings (boundary conditions are deduced from ( |2.8D and matching condition from 

d,iz^A>)U=R, = (2.9) 
a,(z'A<)|,=K = (2.10) 
A<|,=^, = A>|,=y (2.11) 

9,(A<-A>)U=,, = ^ (2.12) 
Solutions of (^^) are linear combination of Bessel Functions: 

A{z',R,R')J^d-i)ipz) + B{z',R,R')N^a-i)ipz) (2-13) 

Applying the conditions (p. 9|)- ( [2^121 ) to this solution leads to an equation which determines 
the KK mass spectrum: 

pR'MpR') + {l-u)J,+i{pR') ^ pRMpR) + {l-u)J,+i{pR) 
pR'N^{pR') + {l-v)N^+i{pR') pRN^{pR) + {l-u)N^+i{pR) ^" ' 

with 1/ = d — 1. Finding the exact solution of ( p.l4| ) is not trivial. To consider only a simple 
case we assume that 4-dim. mass of the scalar field \p\ <C = i.e. pR <C 1. Regarding the 
Standard Model, this can be applied to a confined Higgs when the scale of compactification 
is much higher than Higgs mass or to a light axion like scalar or to graviton with a small 
mass due to radiative corrections. We keep only lowest powers of pR in the expansion of J,^. 
For solving hierarchy problem R' ^ R. Using the asymptotic expansion of Bessel functions, 
(|2.14|) reduces to: 



. — ^ = (2.15) 

1/2-'^ sin(z/7r)r(-z/) ^ ' 
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where ry = pR and rj' = pR'. For > 1 the solutions of ( p.l5|) lead to the following mass 
spectrum: 

ipQ = — iT (2-16) 

m4 ^xyj2v{v - 1) « ^ , T msVS// (2.17) 

v2 

Ip^I fie^'^^imr + + ^) ^^^Se n. (2.18) 

Continuity properties of Ji, guarantees that ( |2.15| ) is also valid when — > or equivalently 
— > 1. In this case r] = or pq = 0. For pR <C 1 the mass difference between KK- modes pn is 
Ap (X 1/ R' <^ 1/ R = fi. Due to special properties of Bessel Function with integer index, the 
zero mode of massless fields is protected from decay even when the spectrum of KK-modes for 
massive particles begins roughly from zero. Tunneling probability depends on 5-dim. mass 
of the field and on warping scale /u. If fi is large, the probability of zero-mode decay to the 
bulk can be small. However, except for very light particles a /i as large as ITeV provides an 
enough large width (F > W~^^eV) for decay to the bulk during propagation in the terrestrial 
atmosphere if m-s is in the mass range of SM particles. 

To see the effect of mass on the coupling we can investigate the mass dependence of the 
propagator on the branes. Using ( |2.8D and ( |2.9D -( ^rT^ ), one can determine the integration 
coefficients A{z' , R, R') and B[z' , R, R') in ( 2.13D and right and left propagators: 



ttR^ I AoJu+i{pz') - AiN^+iipz') i A2Jy+i{pz) - A^N^+iipz) 

2;'^(AoA3 - A1A2) 



ttR^I l 



A2Ju+iipz') - AsN^+iipz') AoJ^+i{pz) - AiN^+i{pz)) 



^>(^'^') = ^ '2(A A ^ A A ^ (2-20) 

2;'^(AoA3 - A1A2) 

Ao = pR'N.ipR') + (1 - iy)N,+iipR') (2.21) 

Ai = pR'MpR') + (1 - i^)J^+i{pR') (2.22) 

A2 = pRN^ipR) + (1 - u)N^+i{pR) (2.23) 

A3 = pRMpR) + (1 - u)Ju+i{pR) (2.24) 

Restricting these equations to the branes gives the 2-point propagators: 

1 f j^iJpi^-x')P'^i^oJu+iipR) - AiN^+i{pR)) 



Aix, R,x',R) = / rfpVP^^-"^ ^ v^u^.+iv^^v (2.25) 

^ ^ {2tt)^J (A0A3-A1A2) ^ ' 

f ™(x-xO?"'(A2 J.+i(pi?') - A3iV.+l(pi?')) 



^ ^ J ^/e-^^-^ - - ^ (AoA:-A,A2) ^'-''^ 



The term R'"^ /R? in ( 2.26D refiects the difference between metric on the branes. Consistently, 



the roots of dominator in ( 2.25| ) and ( |2.26| ), are the same as ( 2.14] ) and correspond to KK- 



modes. Near each mass mode the propagators can be written as a Yukawa propagator with 
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complex mass and a coupling equal to the residue of the integrand. Applying this procedure 
to (|23^ ) one finds that on the brane at R' for z/ > 1: 




£0 



2i 



iv-1) 



(2.27) 



n 



In our approximation |po/^| oc m^/ ^ <^ 1, |p^//u| > |po/mI and it can be even larger than 1. 
Therefore coupling to KK-modes is larger than to zero-mode. The probability of production 
of zero-mode with respect to n''^ KK-mode is: 



For KK-modes with \p'^\ ~ fi the branching ratio V < 1 and therefore the probability of 
production of these modes is larger than zero-mode. 

In conclusion, radiative corrections that can induce an effective mass for scalar/spin-2 
fields weakens their confinement on the brane. In warped 2-brane models even when = 
the mass difference between the zero mode and massive KK-modes on the brane with smallest 
warp factor is small and they can be abundantly produced in high energy interactions. 

2.2 Fermions and Gauge Fields 

Fermions can not be localized to the TeV scale brane (i.e. brane at R') gravitationally but a 
chiral symmetry breaking can confine them [14] . Their escape to the bulk when their > 
has been studied in detail for one brane models in [11] and [29]. The detail of the formalism 
is very similar to the case of a scalar field except that the mass eigenstate equation includes 
a chiral mass term due to interaction of fermion modes with a bulk scalar responsible for the 
symmetry breaking. The width of the zero mode depends on the coupling between fermions 
and the inferred scalar field and without detail knowledge of underlying particle physics it 
is difficult to assess the probability of decay to the bulk. For 2-brane models the situation 
should not be very different and general conclusions of Sec.^]| must be applicable. 

At present the particle physics models don't fix the scale of the symmetry breaking. For 
not creating a new hierarchy however it can not be much larger than compactification scale [18] 
or fundamental scale of gravity M5. Therefore, at energy scale of UHECRs interaction in the 
atmosphere, not only it is possible to produce KK-modes, it is very probable that at such 
energies the restoration of chiral symmetry completely removes the confinement of fermions 
and open the extra-dimension even to fermionic zero modes (presumably SM matter). 

As for gauge bosons, the most successful scenario for their confinement on the branes is 
based on adding an induced kinetic term to the action of bulk gauge fields on the brane. It 
appears due to the interaction of these fields with confined charged scalar or fermions on the 
branes [17]. Other suggestions are mostly equivalent to this scenario [30]. Once the charged 



Po 



2(u-l) , 2 



(2.28) 
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fields become able to escape to the bulk, they drag then' interaction vertex with gauge fields to 
the bulk and release them from confinement. It has been shown [31] that the coupling of gauge 
boson KK-modes to fermions on the brane is stronger than coupling of their corresponding 
zero- mode (similar to the self coupling of scalars ( 2.271 )). 



The general conclusion of this section is that regarding: 

- Very high CM energy of interaction of most energetic Cosmic Rays in the atmosphere 

which is much higher than natural confinement scale of Standard Model particles on the 
brane and natural fundamental scale of gravity to solve the hierarchy problem; 

- The fact that confinement of SM fields is not intrinsic but the result of either a broken 

symmetry (for fermions) or interactions (for bosons); 

- That the particle physics in 5-dim can not be completely massless and at least part of 

the particle spectrum must acquire mass as it is the case in observable 4-dim Universe. 
In fact as radion is in fact the scalar component of 5-dim. metric perturbations [9], 
it couples to all bulk fields and radiatively induces a small mass term. Consequently, 
zero-modes (presumably SM particles) are not stable and in a finite time decay to the 
bulk unless another phenomenon like symmetry breaking prevent it; 

the phase space of production of bulk modes at high energy tail of UHECRs spectra seems to 
be higher than confined modes. 

Until now more than 100 coherent air showers have been observed with Ecm ^ SOOTeV 
(assuming interaction with nucleons in the atmosphere), 17 between them have energies more 
than 450rey and one has a CM energy close to lOOOTeF [32]. Assuming that UHECRs 
interaction with Eqm ^ 300Tey produces bulk modes abundantly, the mere observation of 
coherent showers up to energies close to lO'^^el/ constraints the parameter space of brane 
models. 

As the assessment of cross-sections and other details are model dependent, in the rest of 
this work we simplify the problem of constraining brane models and consider only the classical 
propagation of the bulk modes. This method has been already used by other authors to study 
some of cosmological consequences of brane models [33] [10]. 

3. Propagation 

The geodesic path of particles in the bulk has been already studied in a number of previous 
works [33] [10] [11] [34]. However, most of them are concerned with the possible acausality 
of paths for observer on a brane and their purpose is to see if it can solve the cosmological 
horizon problem in the early universe. Here we are concerned with the present evolution of 
the Universe and simplifying assumptions will be based on its present very slowly changing 
state. 

The metric of the 5-dim brane models can be written as the following: 

ds^ = n^{t,y)dt^ - a^{t,y)5ijdx'dx^ - b'^{t,y)dy^. (3.1) 
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For a static bulk b^{t,y) is constant and we can normalize coordinates such that 6 = 1. The 
geodesic of a particle is defined by: 



-u"u" H u^u^ + 

n n 


i is: 





(3.2) 


du^ 2a' ,• 4 
dr a 


2d ,• n 
a 





(3.3) 


diA 

tin / 

— h nn u u — 

dr 


aau^u-'6ij = 





(3.4) 


dt j dx* 
dr' dr ' 


4 dy 
" =d^- 




(3.5) 



r is the proper time parameter along the particle world line. It is easy to see that (^) is 
integrable and: 

u' = ^ (3.6) 

where 6^ is an integration constant. As we are only interested in the minimum delay in the 
arrival of particles due to the propagation in the extra-dimensions, we put 9^ = 0. Later we 
try to estimate qualitatively the effect of a non-zero 6^. Even after this simplification the 
system of equations ( ^.2D -(^^) is highly non-linear and coupled. In the following we calculate 
an analytical solution for the case h/n ~ 0. This approximation is justified when we are 
interested in the propagation of particles in an extremely short period of time with respect 
to the expansion rate of the bulk or the brane. In fact from the solution of the Einstein 
equations [35], n{t,y) can be normalized such that: 

n{t,v) = % (3.7) 

ao{t) 

where ao{t) = a{t,y = 0) assuming that one of the branes is at y = 0. At present, both 
d{t,y) and do{t) are very slowly varying quantities. Therefore n{t,y) ~ n{to,y) where to is 
the present time, and h{t, y) ~ 0. 
Under this approximation: 



and: 



dr dr 
The solution of (|3.9D and ( 3.10| ) is straightforward: 



dr 



du^ 2u^ dn , , 

— + ——=0 3.9 
dr n dr 

du^ n(u^)'^ dn 

+ ^7^— =0 (3-10) 



= en-\ = ± ( r/ - -2 ) (3.11) 
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The parameters 9 and rj are integration constants and must be determined from the initial 
conditions. The it sign defines the direction of the propagation. In the rest of this letter we 
neglect the direction and consider only the absolute value of u^. 
For a particle leaving the visible brane placed at y = yb at time t = to- 

02 rl Massive particles, 

= u'^{to,ybW{to,yb), V = -tt; n - i^'^ito, yb)f = I 

n [to,yb) Lo Massless particles. 

(3.12) 

After eliminating the proper time from and u^, one obtains the equation of motion in the 
bulk (For simplicity we assume that Dy/dr ~ dy/dr and Dt/dr ~ dt/dr): 



a\t,y) 


= n^{t,y)al{t) = 


A{t) 


= ao\t)-C{t), 


m 


= -/'f,o«o^(i), 


C{t) 


2d2(t) 








/2k2 


1^ 





4a2(t,y) 



dy n^(t,y) 6^ , , 

JL = e 3.13 

dt e \l n^{t,y) ^ ' 

1 Massive particles, 

(3.14) 

Massless particles. 

Our approximations are valid only when dy/dt is real. This put limits on the testable part 
of the parameter space of the models (see below). 

Einstein equations give the solution for a{t,y) and n{t,y) [35]. For a flat visible brane: 

a^(t, y) = A{t) cosh(/iy) + B{t) sinh(/iy) + C(t), (3.15) 

A{t) cosh(/iy) + B{t) sinh(;uy) + C{t)\ 

^, (3.16) 

(3.17) 
(3.18) 

(3.19) 
(3.20) 

For any density />, p' = p/Ajis, J^rs = 'ip/^?- The densities /j^^ and ps are effective total 
energy density of the brane at y = and the bulk respectively. We consider only AdS bulk 
models with pB < 0. The constant = Stt/M^ is the gravitational coupling in the 5-dim. 
space-time. The model dependent details like how p'^^ and ps are related to the field contents 
in the bulk and on the brane and how they evolve are irrelevant for us as long as we assume 
a quasi-static model. The solution ( p. 15 ) is valid both for one brane and multi-brane models. 



The only difference between them is in the application of Israel junction conditions [15] [16]. 
Equation ( p.l4| ) is non-linear and its integration non-trivial. We use again the quasi-static 
properties of the present Universe and its low energy density to simplify the integration. 
Matter density on the branes at late time is much smaller than the brane tension or induced 
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tension by scalar fields [38]. Therefore, it is not unreasonable to neglect time dependence of 
densities in ( |3.17 )-( 3.19| ) and to consider only cosmological constant type energy-momentum 
densities. This simplification is even more justified in our case where we have to deal only 
with very short duration of the propagation in the extra-dimensions. This approximation and 
(|3T7|)-(|3Tg|) lead to: 



c 

A 
B 



2ao{t)ho{t) -C{t), 
-2p'i,^ao{t)ao{t). 



After changing variable y to z = e 



V 



2z 

[l- 



and using (|3?7|): 

V{t,z), 



P'bo 



C(t), 2 2C(t) ,^ , 
"o "o 



C{t). 



(3.21) 

(3.22) 
(3.23) 



(3.24) 
(3.25) 

(3.26) 



f =,yP(.-^P). 

If an ejected particle to the bulk comes back to the brane, must go to zero at some point 
in the bulk before the particle arrives to the bulk horizon (if it is present). The roots of 
( |3.26| ) correspond to these turning points and determine the propagation time in the bulk. 
In the next simplifying step we use again the fact that the typical propagation time we are 
interested in is very much shorter than the age of the Universe and therefore A, B, C and 
A,B,C during propagation are roughly constant, the right hand side of (|3.26|) depends only 
on z and is easily integrable: 

/■"^"'"p 2dz 

^y^stop 



At 



propag 



to) 



20 



^ V{z--^V) 



(3.27) 



In (3.27), to is the initial time of propagation in the extra-dimension and tstop is the time 
when the particle's velocity changes its direction, i.e. when dz/dt = 0. The integral in (3.27) 



is related to the elliptical integrals of the first type J-{ijJ, v) where to and v are analytical 
functions of the denominator roots in ( p. 27 ) and zq [39]. Note that Zgtop corresponds to the 
closest root to zq. 



4. Test of Brane Models 

In this section we apply the formalism discussed in the previous section to most popular brane 
models and determine the propagation time of high energy particles in the fifth dimension. 
Note that the calculation of propagation time in these models under our approximations is 
valid only for durations very smaller than the age of the Universe and if in the following 
figures in part of the parameter space the propagation time can be larger, this part of the 
figure should not be considered. 
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4.1 2-Brane Models 



It has been shown in [15] [16] that by 
imposing constraints on the visible brane to 
obtain the observed value of cosmological 
constant A and Newton coupling constant 
G and to solve the hierarchy problem, all 
parameters of this class of models i.e. po, 
PL and h.Rs, can be determined as a func- 
tion of fiL where L is the distance between 
two branes. It is not however possible to 
find an exact analytical form for the solu- 
tions. Moreover, the analytical solution in 
[16] has been obtained for a special setup 
which decouples hidden and visible branes. 
Here we free some of the constraints, first to 
be able to find analytical solutions, and sec- 
ond to extend this study to a larger number 
of models. 

4.1.1 Geodesies in Static RS Models 

In the original RS model with static metric: 




20 



40 
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80 



1 00 
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Figure 1: Propagation time for relativistic particles 
with ul{tQ)/N = 10^ (full line) and <(to)/A^ = 
1.2 (dash line) in RS model. Red, magenta and 
light green curves correspond to A/5 = lO^^eT^, 
Ms = lO^^eV and M5 = lO^^eV (or 11 ~ IQ-^eF, 
/i ^ lO^^^e^ and /i ^ IQ^^eV for fine-tuned model) 
respectively. The dark green line shows the time co- 
herence precision of present Air Shower detectors. 



dy\ 



(4.1) 



the main constraint on the model is the cancellation of the cosmological constant on the visible 
brane which leads to the equal and opposite sign tensions pQ = —pi = Ars- The solution of 
the hierarchy problem limits the range of the parameter pL. It has been shown [10] that in the 
fine-tuned RS model photons leave the brane and never return. Assuming the visibility of the 
all dimensions of the space-time at very high energies, in this model we could never observe 
the ultra high energy particles and therefore it is automatically ruled out. Nonetheless, to 
test the formalism of the previous section we apply it to this model. 

For a very small cosmological constant (as it is assumed in the process of fine-tuning) C{to) ~ 
and 'D{tQ, z) ~ 1. For massive particles, the denominator in ( |3.27| ) has only one root: zi = 1/9 
and: 



At, 



propag 




(4.2) 



Fig.|l] shows At function of pL and M5 for massive relativistic particles. With 

present air shower detectors time resolution of order lO^^sec, only when M5 > lO^^eV, the 
model is compatible with the observed time coherence of the UHE showers. For fine-tuned 
RS model p « G/k'^ [3] i.e. p = GM^ ~ IQ-'^eV for M5 ~ lO^^eF [3]. Due to smallness of 
p and consequently lightness of KK-modes for SM particles even this model with large M5 
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has already been ruled out [31] unless a conserved quantum number prevents the production 
of KK-modes [21]. 

For massless particles: 

z{t) -zo = fi{t - tof (4.3) 

In ( [4. 3D , z{t) is monotonically increasing and there is no stopping point. With our approxi- 
mations there is no horizon in the bulk because a{t, y) is roughly constant. Therefore ( [4. 3D 
means that massless particles simply continue their path to the hidden brane and their faith 
depends on what happen to them there. At very high CM energy of UHECR interactions if 
charged particles can escape to the bulk photons are also dragged to the bulk and never come 
back. 



4.1.2 General Solution 



Numerical solution of constrained 2-brane models in [15] [16] shows that for fiL > 5 the 
tension on both branes is positive and very close to A/jg. We can use constraints on the 
Cosmological Constant and hierarchy to find p'q and p'j^ = p'^. We redefine them as p'^ = l+Apg 
and = 1 + Ap^. To solve hierarchy problem (See equations 29-31 in [16]): 



Ml 



= !!i - ^Aq^^ ~ cos^^(/^-^)) + sinh(/iL) 
rig (1 — cosh(pL)) -|- sinh(pL) 



<. 1 



This leads to: 



Apl 



iV2 ^1 _ g-Mi + Ap'^(l - cosh(/iL)^ ^ 



1 — cosh(/iL) 



1 — cosh(pL) 



For a very small A^^ and Ap^ < 1, the first term in (|4.5|) is 0{N'^) and: 



Ap'o 



1 - e-^^ 



2e 



1 — cosh(/iL) 1 — cosh(/iL) 
Using a\/a?j^ = where H is the Hubble Constant on the visible brane [16]: 



2A2 sinh(pL) 



(1 — cosh(/iL)) 



2/?2 



+ e 



1 ± 



2H^ 2H^ 
((1 - cosh(pL))^^ + e-^'^ - 1)2 + A2 sinh(;uL)(^^ + 2e-^'^] 
p^ p^ 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



In {'L7) the solution with plus sign gives Ap^ ~ —2 which deviates from our first assumption 
|Ap^| < 1 and leads to a negative tension on the visible brane like static RS model. The 
solution with negative sign is: 

Ap'l « 2e-''^ (4.8) 
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and both branes have positive tension close to A^s- 

When the matter densities on the branes and in the bulk are negligible [16], Oq/oq = a\/a?^ 
H'^ and: 



C{t) 



= c' 



2H' 



(4.9) 



It is easy to see that V and a^(t, y) have the same roots. Models with a horizon i.e. a point yh 
such that a'^{t,yh) = are pathological (because no particles/brane behind it is observable). 
The condition to have no real root i.e no horizon in the bulk is: 



Q/2 Q/2 

-2 + —^Ap'o + C'ii- — 
For massive particles, the denominator of the integrand in ( p. 27 ) can have two roots: 

c'-e^± V(C' - 02)2 + (2 + ap'q + c'){Ap'q + c) 



(4.10) 



z± 



(4.11) 



The model is consistent only if P(z — > in the range of integration. Therefore: 



(4.12) 



The matter on the brane is confined only if 2:_|_ > 1 and — > zq when ^ 0. To first order 
in e~'*^ and A^^ ^^ig leads to the following relation between parameters of the model: 



-C' = Ap', + -iN' + Ap',) 
zq 



(4.13) 



This condition is not an addition to the model described in [16]. It is in fact the result of 



solutions (4.6) and ( [4.^ ) under the approximations considered here. It is not evident whether 
such a constraint appear in the full theory. 

For pL > 5 the right hand side of (4.13) is positive. Therefore —C oc can not be zero. 
This relation between a small but non-zero value of the Hubble Constant or equivalently Cos- 
mological Constant on the visible brane and the smallness of N'^ and p which is related to the 
strength of the induced gravitational coupling on the brane, confirms the same observations 
in [16] for an analytical solution of 2-brane models with some approximations and in [37] 
for the exact solution of some special models. 
Finally the propagation time in the bulk is given by: 



At 



propag 



a 



M8|Ap'o + C'|)4 
2 arctan 



-^(«,Q). 



iq(z. 



Zo) 



p{zo - z+) 



Q 



'2 + 



pq 



C'iC - 02) + 4(2 + Ap'q + C'){Ap'q + C) 



(Ap'. + C 



/^2 



(4.14) 
(4.15) 
(4.16) 
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c 



C 



2 + Apf, + C 



C 



Ap'^+C 



Figj2| shows the propagation time for models which satisfy simultaneously (|47 



(4.17) 
and 





Figure 2: Left: Propagation time for relativistic particles in 2-brane model. Description of the curves 
is the same as FigJ^. Right: Parameter y.[eV) as a function of fiL. It is roughly independent of M5 
(See the text). 



( [4.13| ). In equation ( |4.13 ) up to first order, C' depends only on fiL and thus in (4.9) the 
value of /X is independent of M5. Only models with large /uL ^150 are not ruled out. This 
is due to ( 4.10| ) and smallness of the observed Hubble Constant H. The same conditions 
make p > leV which is much higher than the for the fine-tuned RS model. With fi < 
rriradion [6] [9] these class of brane models are also consistent with constraint on the fifth- 
force measurements [41]. From ( ^.17 ) the corresponding life-time for SM fields with ms > 
is shorter than W~^^sec, much shorter than propagation time in the atmosphere and also 
much shorter than propagation time in the bulk. This justifies the classical treatment of 
propagation. 

We have also applied the formalism described here to the fine tuned model of [16]. In 
this model the equation of state on the hidden brane is fine-tuned to neutralize its effect on 
the visible brane. This results to a unique definition of G and the only free parameter in the 
model is M5. Although this model has been obtained just by phenomenological arguments, a 
field theory model suggested by Arkani-Hamed et al. [36] to solve the Cosmological Constant 
problem has the same form for G if fxL ~ k^(/)(0) where </>(0) is the vev of radion on the visible 
brane (The Arkani-Hamed et al. model has only one brane but it includes a horizon in the 
bulk which limits the accessible size of the extra-dimension and makes it similar to a two 
brane model). 
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Figj3| shows the propagation time and ^ clS cl function of M5. Unfortunately despite 
physical interest of this model it is only compatible with very high M5 > lO^^el/ unless the 
decay to the bulk is prevented up to such high energies. The corresponding // however is 
consistent with present constraint on the Fifth force [41] and the life-time of zero-modes of 
massive 5-dim. fields for such high M5 is enough short to permit particles to decay to the bulk 
during their propagation in the atmosphere. The value for M5 is some 6 orders of magnitude 
larger than WTeV, presumably the Electroweak interaction scale and it would be a matter 
of speculation to consider this model as having no hierarchy problem. Another problem in 
testing these models with UHECRs is that as the natural scale of gravity M5 is very high, it is 
possible that the symmetry breaking scale which is necessary for the localization of fermions 
(and indirectly gauge bosons) is also much higher than CM energy of UHECRs interaction. 
In this case only very weakly interacting particles like gravitons can decay to the bulk. As 
the total production cross-section for them can be tiny, the number of observed UHECRs 
event can be not enough to constrain such models. We have also tested the general 2-brane 





bg M5 log M5 

Figure 3: Left: Propagation time for relativistic particles in the fine-tuned 2-brane model of [16]. 
Right: Parameter as a function of M5. 



models without taking into account ( 4.13 ). Roughly speaking, it is equivalent to having a 
comparable matter density and tension on the hidden brane. The value of iJ, becomes a free 
parameter. The result is shown in Fig.^ for 3 different values of fi. Models with /j, > lO^eT^ 
and fiL < 70 are compatible with the present observation of UHECRs. The lower limit for /i 
from this test is higher than the constraint obtained from Fifth force experiments [41]. 

4.2 One-Brane Models 

The solution of Einstein equations for symmetric one-brane models is the same as two-brane 
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Figure 4: Propagation time for relativistic particles in 2-brane models with /i as a free parameter. 
Top left: ^ = IQ-'^eV; Top right: ^ = lO^'ey; Bottom: ^ = lO'^eV. Description of the curves is the 
same as Fig.|l|. 



ones [35]. Due to existence of only one boundary however the bulk and brane tensions are 
not related. In addition, the cosmological evolution on the brane: 

^ = ^PB + ^{pt + Pm{t)f + ^. (4.18) 

includes an arbitrary function C (t) which is related to the bulk tension and matter [40] . Here 
we test two popular models studied in [35] and [40]. 

In the first model [35] C{t) =0 and the brane tension is fine-tuned to cancel the effect of 
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quadratic term at late times: 



and: 



SvrG 



6 



(4.19) 



(4.20) 



The only free parameter in the model is M5. In the second model [40] C{t) (or equivalently 
Tg component of the energy-momentum tensor) is adjusted such that the conventional evo- 
lution equation be obtained. At late times when the brane tension is much larger than time 
dependent matter terms these two models are roughly the same. 

Equations ( 4.19|) and ( 4.20 ) determine pb and fi. It is easy to see that ph = Ajis- The 
definition of C and roots are the same with Apo = ^Pb = and zq = 1. Fig.^ shows the 
propagation time for these models. 
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Figure 5: Left: Propagation time for one-brane models of [35]. Right: Parameter /j, as a function of 
pL. 



4.2.1 Effect of 61* / 

When 9^ ^ equations (| 



and ( 3.10 ) can be written as: 



du° 2n0 dn hO'eHr 



+ 



+ 



u 



dr n dr 



dr 



At least formally Eq.( |42TD can be solved analytically 

_ e'e^,,g{t,y) 



— + 



(4.21) 
(4.22) 

(4.23) 
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dr 



(4.24) 



We don't need to solve ( [4.22| ) directly. Knowing and we can use the definition of velocity 
vector to determine n'^: 

n\u°f - a^VSij - {u^ = e (4.25) 

The definition of e is the same as in (|3.14 ). After elimination of dr we obtain the formal 
description of the equation of motion in the bulk: 



dt 



+ 



We can use ( |3.6| ) to determine dr (As before for simplicity we assume that Dx"^ /dr 

dr 

and: 



a^dijO^dx' 



e'e^i,Q{t{T),y{T)) 



-SijO^dx' 



(4.26) 

dxydr): 
(4.27) 

(4.28) 



In ( 4.28| ) d/a is independent of x*. The rest of right hand side of ( [4.28] ) i.e. J SijO^dx'^ 
is the projection of the particles world line on the brane. The value of Hubble constant 
d{t,y)/a{t,y) ~ H'^ (from ( ^TploD and (|]2jg)) and thus 9'9^5ijg is very smah when 
the projection distance traversed by the particle is small with respect to the Hubble radius. 
Therefore we presume that conclusions of the previous section will not be extremely modified 
when full propagation is considered. 



5. Conclusion 



The calculation in this work is mainly based on two assumptions: 

- At interaction energy scale of most energetic Cosmic Rays the physics is high dimensional 

either because all symmetry based confinements are no longer at work or because there 
is a large cross-section and/or phase space volume which permits the production of 
massive KK-modes. 

- In the time scale of the propagation of a particle in the extra-dimension, the bulk and the 

branes are quasi-static. 

If these assumptions are valid the time coherence of Ultra High Energy Air Showers rules out 
a large part of the parameter space for a number of brane models unless some micro-physics 
phenomena confine particles to the brane at energies much higher than Electroweak scale. 
This makes a new hierarchy inconsistent with the spirit of brane models. 
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For most 2-brane models the acceptable range of ;U is /i > leV except for original RS 
model which needs /x > 10~'^eV. This lower limit is much lower than what can be obtained 
from non-observation of KK-mode production in accelerators [31]. However, at present ac- 
celerator energies it is always arguable that presence of some quantum conservations prevent 
the production of KK-modes. Presence of such conservations at the CM energy of UHECR 
interaction seems much less natural and constraints on M5 more arc robust. 

The upper limit of L is also a universal value for models with different range of //L: 
L < lO'^eV"^ ~ 10~^cm. It is much smaller than the upper limits obtained from gravity 
experiments [42] [41]. Again for static RS model the upper limit is L < 10^ cm, much larger 
than one obtained for other models. One brane models with interesting range of M5 are ruled 
out. 

This study has an additional interesting conclusion: The close relation between a very 
small but non-zero Cosmological Constant and the smallness of the Newton coupling constant 
(i.e. the hierarchy problem). In fact without the fine-tuning of these apparently independent 
physical quantities, the brane models are not consistent. 

References 

[1] Kaluza Th., Sitz. Preuss. Akad. Kl, 966, 1921, Klein O., Z. Phys. 37, (1926) 895. 

[2] Arkani-Hamed N., Dimopoulos S. & Dvali G., Phys. Lett. B 429, (1998) 263 hep-ph/9807344; I. 
Antoniadis, et al. , Phys. Lett. B 436, (1998) 257 hep-ph/9804398. 

[3] Randall L. & Sundrum R., Phys. Rev. Lett. 83, (1999) 3370 hep-ph/9905221, Phys. Rev. Lett. 
83, (1999) 4690 hep-ph/9906064. 

[4] Rubakov V.A. & Shaposhinkov M.E., Phys. Lett. B 125, (1983) 139, Akama K., Lect. Notes 
Phys. 176, (1982) 267. 

[5] Horava R & Witten E., Nucl. Phys. B 460, (1996) 506 hep-th/9510209, Nucl. Phys. B 475, 
(1996) 94 hep-th/9603142. 

[6] Goldberger W.D., Wise M.B., Phys. Rev. Lett. 83, (1999) 4922 hep-ph/9907447. 

[7] Binetruy R, Cline J.M. & Grojean C., Phys. Lett. B 489, (2000) 403 hep-th/0007029. 

[8] Kobayashi Sh., Koyama K. & Soda J., Phys. Lett. B 501, (2001) 157 hep-th/0009160. 

[9] Csaki C, Graesser M., Randall L. Terning J., Phys. Rev. D 62, (2000) 045015 hep-ph/9911406, 
Charmousis Ch., Gregory R. & Rubakov V.A., Phys. Rev. D 62, (2000) 067505 hep-th/9912160. 

[10] Gregory R., Rubakov V.A. & Sibiryakov S.M., Phys. Rev. Lett. 84, (2000) 5928 

hep-th/0002072, Miick W., Viswanathan K.S. & Volovich I.V., Phys. Rev. D 62, (2000) 105019 
hep-th/0002132, Gregory R., Rubakov V.A. & Sibiryakov S.M., Class. Quant. Grav. 17, (2000) 
4437 hep-th/0003109. 

[11] Dubovsky S.L., Rubakov V.A., Tinyakov RG., Phys. Rev. D 62, (2000) 105011 hep-th/0006046, 
J. High Ener. Phys. 0008, (2000) 041 hep-ph/0007179. 



-20- 



[12] Jackiw R. & Rebbi C. Phys. Rev. D 13, (1976) 3398, Jackiw R. & Rossi R, Nud. Phys. B 190, 
(1981) 681, t'Hooft G., Phys. Rev. D 14, (1976) 3432, Visser M., Phys. Lett. B 159, (1985) 22 
hep-th/9910093. 

131 Gherghetta T., Roessl E. & Shaposhinkov M., Phys. Lett. B 491, (2000) 353 hep-th/0006251. 

141 Bajc B. & Gabadadze G.,Phys. Lett. B 474, (2000) 282 hep-th/9912232. 

151 Kanti P., Olive K.A. & Pospelov M., Phys. Rev. D 62, (2000) 126004 hep-ph/0005146. 

161 Ziaeepour H., "Two-Brane Models and BBN" hep-ph/0010180. 

17] Dvali G., Gabadadze G. & Shifaian M., Phys. Lett. B 497, (2001) 271 hep-th/0010071, 
Dubovsky S.L., Rubakov V.A., Int. J. Mod. Phys. A 16, (2001) 4331 hep-th/0105243. 

18] Carena M., Tait T.M.P., & Wagner C.E.M., hep-ph/0207056. 

191 Dvali G., Gabadadze G. & Porrati M., Phys. Lett. B 484, (2000) 112 hep-th/0002190, Phys. 
Lett. B 484, (2000) 129 hep-th/0003054. 

201 Dubovsky S.L., Rubakov V.A. & Tinyakov P.G., Phys. Rev. Lett. 85, (2000) 3769 
hep-ph/0006046, J. High Ener. Phys. 0008, (2000) 041 hep-ph/0007179. 

21] Appclquist T., Cheng H.C. & Dobrcscu B.A., Phys. Rev. D 64, (2001) 035002 hep-ph/0012100. 

Servant G. & Tait T.M.P., hcp-ph/0206071 

22] Dimopoulos S. & Giudicc G.F., Phys. Lett. B 379, (1996) 105, Long J.C., Chan H.W. & Price 
J.C., Nud. Phys. B 539, (1999) 23 hcp-ph/9805217. 

23] Arkani-Hamed N., Dimopoulos S. & Dvah G, PLB 429, (1998) 263 hep-ph/9807344, Fairbairn 
M. & Griffiths L.M., J. High Ener. Phys. 0202, (2002) 024 hep-ph/0111435. 

24] Giudice G.F., Rattazzi R., & Wells J.D., Nud. Phys. B 544, (1999) 3 hep-ph/9811291. Hall L.J. 
& Smith D., Phys. Rev. D 60, (1999) 085008 hep-ph/9904267. 

25] Davoudiasl H., Hewett J.L. & Rizzo T.G. Phys. Rev. D 63, (2001) 075004 hep-ph/0006041. 

26] Fairbairm M. & Griffiths L.M. J. High Ener. Phys. 0202, (2002) 024 hep-ph/0111435. 

27] Giddings S.B., Katz E. & Randall L., /. High Ener. Phys. 0003, (2000) 023 hep-th/0002091. 

28] Rubakov V.A., Phys. Usp. 44, (2001) 871 hcp-ph/0104152. 

29] Ringcval Ch. Peter P. & Uzan .J.P., Phys. Rev. D 65, (2002) 044016 hep-th/0109194. 

30] Akhmedov E.Kh, Phys. Lett. B 521, (2001) 79 hep-th/0107223, Tachibana M., hcp-th/0108164. 

31] Davoudiasl H., Hewett J.L. & Rizzo T.G., Phys. Lett. B 473, (2000) 43 hep-ph/9911262. 

32] Bird D., et al. , Astropart. Phys. 441, (1995) 144, Takeda M., et al. , Proc. 27th Int. Conf. on 

Cosmic Rays, Hamburg, 2 (2001) 675. 

33] Chung D.J.H. & Frccsc K., Phys. Rev. D 62, (2000) 103502 hcp-ph/9910235. 

34] Davis A.Ch., Rhodes Ch. & Vernon I., J. High Ener. Phys. 0111, (2001) 015 hep-ph/0107250. 

35] Binetury P., DefFayet C. & Langlois D., Nud. Phys. B 565, (2000) 269 hep-th/9905012, Phys. 
Lett. B 477, (2000) 285 hep-th/9910219. 



- 21 - 



[36] Arkani-Hamed N, Dimopoulos S., Kaloper N. & Sundrum R. Phys. Lett. B 480, (2000) 193 
hep-th/0001197. 



[37] Flanagan E.E. et al. , Phys. Rev. D 64, (2001) 045007 hep-th/0012129, Flanagan E.E., Tye H. 
& Wasserman I., Phys. Lett. B 522, (2001) 155 hep-th/0110070. 



[38] Brax Ph. & Davis A.C., hep-th/0104023, Davis S.C. J. High Ener. Phys. 0203, (2002) 058 
hep-ph/0111351. 



[39] Gradshteyn I.S. & Ryzhik I.M. "Table of integrals, series and products" , Academic Press INC. 
1965. 



[40] Kanti P., Kogan 1.1., Olive K.A. & Pospelov M., Phys. Lett. B 468, (1999) 31 hep-pli/9909481, 
Phys. Rev. D 61, (2000) 106004 hep-ph/9912266; Csaki C. et al. , Phys. Rev. D 62, (2000) 
045015 hep-ph/9911406. 



[41] Will CM. Living Rev.Rel. 4, (2001) 4 gr-qc/0103036, Bae S. & Seok Lee H. Phys. Lett. B 506, 
(2001) 147 hep-ph/0011275, Perivolaropoulos L. & Sourdis C. hep-ph/0204155. 



[42] Long J.C., Churnsidc A.B.& Price J.C., Proceedings of the Ninth Marcel Grossmann 
Conference (Rome, 2-8 July 2000) hep-ph/0009062. 



- 22 - 



